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Abstract 



We present a detailed classification of random Dirac hamiltonians in two 
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teen classes. We also extend this classification to non-hermitian hamiltonians 
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I. INTRODUCTION AND RESULTS. 



Recently, there is an increasing interest in two dimensional localization problems whose 
behavior differs from generic Anderson localization [1]. To mention but a few examples, 
there are investigations of the quantum Hall transition [2-5], of quasi-particle localization 
in systems with degenerate Fermi surfaces [6] and dirty superconductors [7-10], and studies 
of hopping models on bipartite lattices [11,12]. 

Universality classes of localization/delocalization transitions depend largely on their dis- 
crete symmetries. Thus the Wigner-Dyson classification of random hermitian matrices [13] 
[14] plays a significant role. Localization in superconductors led Altland and Zirnbauer 
to significantly extend the Wigner-Dyson classification by incorporating particle-hole and 
chirality symmetries of the matrices [15]. 

Most of the localization problems mentioned above may be formulated as spectral prob- 
lems for Dirac-like hamiltonians in two spacial dimensions, and many of them differ from 
generic Anderson localization in that they exhibit a singular density of states at the critical 
point. 

In this note, we present a somewhat detailed classification of such Dirac hamiltonians in 
two dimensions. This classification potentially differ from that of random matrices because 
in the latter no structure is imposed on the matrices — and matrices differing by unitary 
similarity transformations are treated as equivalent — whereas in the former the hamiltoni- 
ans a priori possess a Dirac form. This may have two opposite effects: either some random 
matrix classes may not be realized by Dirac operators, or Dirac hamiltonians belonging to 
the same random matrix class may not be equivalent if the unitary similarity transforma- 
tion relating them does not preserve the imposed Dirac structure. Although some models 
of random Dirac fermions have already been identified with Altland-Zirnbauer classes, the 
correspondences have not been fully established and one motivation of our classification is 
a more complete dictionary. 

We consider Dirac hamiltonians H = (r x p x + r y p y )/2 + f ■ W + Wo where f are Pauli 
matrices 1 , p x<y = —id XtV , and Wo and W are generalized masses or potentials, and are 
matrices acting on an isopsin sector. Introducing complex coordinates, z = x+iy, z = x—iy, 
and d z = (d x — id y )/2, = (d x +id y )/2, after a unitary transformation one generally obtains 
the following 2x2 block structure: 



Here A z , A- and V± are random matrices depending on the spacial coordinates x, y, and 
belonging to some statistical ensemble. 

As usual, the classes are sets of hamiltonians with specific transformation properties 
under some discrete symmetries. For Dirac hamiltonians (1.1), the simplest symmetries are 
chiral, particle-hole, and time-reversal symmetry, which relate the hamiltonian H to —H, its 



Pauli matrices will be denoted f or a depending on which space they are acting. Our convention 
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transpose H T and its complex conjugate H* respectively. We demand that these transforma- 
tions are implemented by unitary transformations and that their actions on the hamiltonian 
square to one. They also should preserve the form (1.1) of the Dirac hamiltonian. Hence 
we consider the following transformations: 



P sym. : 
C sym. : 
K sym. : 



H = -PHP-\ 



H = e k KH*K-\ C 







-1 J 
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-e c a 


o 





K 


-e k K 






pp ] = 1, p 2 



KK^ = 1, K T = ±K 



(1.2) 
(1.3) 
(1.4) 



where e c = ±1. Type P symmetries are commonly referred to as chirality symmetries, C 
expresses a particle-hole symmetry, and K time-reversal symmetry. For hermitian hamilto- 
nians, since H T = H*, C and K symmetries are identical and we will only talk about C 
symmetry, where e c = +1 will be interpreted as time- reversal symmetry and e c = — 1 will 
be referred to as particle-hole symmetry. 

We found thirteen distinct classes of hermitian Dirac hamiltonians, listed in eqs.(2.16- 
2.25). This classification, which is presented in Section II, is slightly finer than that of 
random hermitian matrices [15]. See Table 1 for a comparison. (The numbering of the 
classes has no special meaning.) The essential difference between these two classifications 
is a doubling of the chiral classes. This arises from differences in the notion of equivalent 
hamiltonians. 





Random matrix 


Time rev. 


Part. -hole 


Chirality 


Sym. 
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inv. 


symmetry 




group 


class 


A=GUE 


no 


no 


no 


U{n) 


class 1 


A m =chiral GUE 


no 


no 
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U(n) 


class 2 


A in =chiral GUE 


no 


no 


yes 


U{n) x U(n) 


class 3+ 


A„=G0E 
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no 


no 


0(n) 


class 3 
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no 


yes 


no 


0(n) 


class 4 + 


Ai=GSE 


yes 


no 


no 


Sp(2n) 


class 4 
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yes 


no 


Sp(2n) 


class 5 


D m =chiral GOE 
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yes 


yes 


0(n) 


class 6 


Ci =chiral GSE 


yes 


yes 


yes 


Sp(2n) 


class 7 


D m =chiral GOE 


yes 


yes 


yes 


0(n) x 0(n) 


class 8 


Ci=chiral GSE 


yes 


yes 


yes 


Sp(2n) x Sp(2n) 


class 9 + 


Di 


yes 


yes 


yes 


U(n) 


class 9 


c„ 


yes 


yes 


yes 


U(n) 



Let us mention a few well-known realizations of the classes of Table 1. Classes (GUE), 
3 + (GOE) and 4 + (GSE) are the usual Wigner-Dyson classes. The U(l) model of class 
was introduced in [5] in connection with the quantum Hall transition. The U(n) case 
of class appeared in [6] for describing systems with degenerate Fermi points. The chiral 
classes 1 (chGUE), 5 (chGOE) and 6 (chGSE) are realized by Dirac operators coupled only 
to random gauge potentials [19]. The U(n) model of class 1 was applied to dirty d-wave 
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superconductors in [7]. The pure U(l) random gauge potential is still not completely solved 
due to some recently recognized non-perturbative effects [16-18]; the non-abelian cases have 
been solved by various methods [7,20,21]. The chiral class 2 (chGUE) is realized by the 
Gade-Wegner hopping models [11]. (See [22] for a recent numerical study.) The class 3 
(D) [23] and class 4_ (C) [24] appeared in the context of dirty superconductors with broken 
time reversal symmetry. Numerical analysis of the sp(2) model of class 4_ (C), the so-called 
spin quantum Hall effect, was performed in [25] and exact results were obtained by mapping 
it to percolation [26]. 

Since each of the ten classes of random matrices has been associated with a sigma-model 
on a symmetric space, an interesting open question concerns how the sigma-models can 
incorporate the finer classification described in this paper. In particular the renormalizable 
effective field theories described in section II can have 1, 2, ... up to 10 couplings whereas 
sigma models on symmetric spaces generally have a single coupling and possibly an additional 
coupling coming from a topological #-term or Wess-Zumino term. 

Non-hermitian random hamiltonians have recently been used in the description of var- 
ious phenomena, see e.g. [27]. In Section III we extend our classification to non-hermitian 
hamiltonians. In addition to type P, C and K symmetries, we may consider a Q symmetry 
relating H to its adjoint: 



Imposing these symmetries selects reality conditions on the potentials A z , A- and V±. Of 
course, the type Q symmetry with e q — 1 and £ = 1 simply means that H is hermitian. 
From any non hermitian Dirac operator one may naturally define a hermitian one, denoted 
H, by doubling Hilbert space on which it acts. (See eq.(3.4).) The latter hamiltonian is 
then a representative of the chiral class we have indexed as class 2 (chGUE). However, as we 
explain, the classification of the non-hermitian Dirac operators H is finer, and more involved, 
than that of their doubled companions 7i. We find a total of 87 universality classes. 

As a byproduct of our analysis we are easily able to classify non-hermitian random 
matrices without a Dirac structure. This leads to 43 classes. 

II. CLASSIFICATION OF HERMITIAN DIRAC HAMILTONIANS. 

We first consider hermitian hamiltonians which requires Vt = V± and A\ = A-. Let 
us define a "minimal class" as a class of Dirac hamiltonians which cannot be simultane- 
ously block diagonalized. For such hamiltonians there exists no fixed unitary matrix S that 
commutes with the hamiltonian, H = S H S -1 , and preserves the Dirac structure, which 
requires S = diag(s,s). The existence of an integral of motion, such as spin, implies such 
an S and the resulting hamiltonian is thus not minimal according to our definition. Rather 
our ensembles apply to each block with fixed quantum numbers of the integrals of motion. 



Only type P and type C symmetries are relevant for hermitian hamiltonians. We first 
need to classify the compatible operators P and C, or preferably the compatible 7 and a. 




(1.5) 



A. Compatible symmetries. 
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It is important to bear in mind that what is meaningful is the group generated by these 
symmetries. For instance, if the hamiltonian possesses both a P and a C symmetry, then it 
automatically has another C-type symmetry C: 

H = e' c C H T C'-\ C = PC, e' c = -e c (2.1) 

For hermitian hamiltonians, since C can be interpreted as a time-reversal (particle-hole) 
symmetry if e c = — 1 (e c = +1), the classes with both a P and C symmetry thus automati- 
cally have chirality, particle-hole and time-reversal symmetry. 

The operators P and C are defined up to dilatations by scalars and up to unitary changes 
of basis, H — > U HW, which preserve the form of the Dirac hamiltonians. This requires 
U = diag(w, u). On 7 and a, this translates into: 

7 -> «7M^ ; a^uau T (2.2) 

with u unitary. The unitarity and the order two constraints on P and C imply: 

77 t = 1 ? 7 2 = l . a(J \ = l a T = ±(J 

These conditions are covariant under the gauge transformations (2.2). 

Let us first only impose a type P symmetry. Modulo (2.2) we can reduce 7 to a diagonal 
matrix with only ±1 on the diagonal. We may thus choose: 

case 1) : 7 = 1 (2.4) 
case 2) : 7 = a z <g) 1 (2.5) 

In the second case, we assumed for simplicity that the numbers of +1 and —1 in 7 are equal, 
but this could be generalized. 

Let us now impose only a type C symmetry. Up to the transformations (2.2), there are 
two (standard) cases [13] depending on the condition a T = ±a: 

case 3) : a — 1 (2.6) 
case 4) : a = ia y <g> 1 (2.7) 

Indeed, assume that a T = a. Then, since a is also unitary, a a* = 1 and its real and 
imaginary parts commute and are both symmetric. They can be simultaneously diagonalized 
by a real orthogonal matrix o, so that a = o5o T with 5 a diagonal unitary matrix. Hence, 
a = uu T , with u = o5 l l 2 unitary, and a is equivalent to the identity modulo (2.2). The 
argument is similar for a T = —o. 

Next we impose simultaneously a type P and a type C symmetries. These symmetries 
have to be compatible in the sense that their actions on hamiltonians should commute. 
For generic hamiltonians this requires that C oc PCP T , or a oc ^a^ T . This condition is 
covariant under transformations (2.2). So we may choose a basis in which 7 is diagonal and 
we restrict ourselves to the two cases (2.4,2.5). We then have two sub-cases corresponding 
to the two possible values ±1 of the proportionality coefficient in the above equation, so 
that a either commutes or anticommutes with 7: 

(j = ±7(X7 T =>• [<7, 7 ] = or {(j,7} = (2.8) 
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If 7 = 1, cr automatically commutes with it and we obtain: 



case 5) : 7 = 1, a = 1 (2.9) 
case 6) : 7 = 1, cr = ia y <g> 1 (2-10) 

If 7 = <t z <S> 1 , we have to consider separately the two possibilities in (2.8). The transforma- 
tions (2.2) have to preserve the form of 7 so that u has to be block diagonal u = diag(«i, w 2 ) 
with Mi )2 unitary. When [cr, 7] = 0, a has also to be block diagonal, a = diag (cri, cr 2 ). As 
above, modulo (2.2) with u = diag (it 1, u 2 ), it can be reduced to a = 1 if a T = cr and 
a = io y <S> 1 if cr T = —a. Thus, we get two possibilities: 

case 7) : 7 = a z <g) 1, cr = 1 2 <S> 1 (2.11) 
case 8) : 7 = cr z <g> 1 2 ® 1, a = 1 2 ® icr y <g> 1 (2.12) 

When {cr, q] = 0, cr has to be block off-diagonal, so that a = ( T *V with s unitary, 



depending whether cr is symmetric or antisymmetric. The gauge transformations (2.2) then 
become s^itiSMj with ui j2 unitary, and any unitary s is gauge equivalent to the identity. 
This gives two cases: 

case 9) : 7 = a z <g> 1, cr = icr^ <g> 1 (2-13) 
case 9') : 7 = o z <S> 1, cr = a x <g> 1 

Cases 9) and 9') turn out to be equivalent because the type C symmetry of one of the two 
cases follows from the product of the type P and the type C symmetries of the other case. 

Finally let us consider more combinations of type P or C symmetries. If we impose two 
symmetries of type P, their product commutes with the hamiltonians and this system thus 
does not correspond to a minimal class. Next consider imposing two compatible symmetries 
of type C with sign e c i and e c2 . If the product e c ie c2 = —1, their product (see eq. (2.1) ) 
makes a type P symmetry. Thus two type C symmetries with opposite e c signs are equivalent 
to a type P and a type C symmetries which we have already classified. If the e c signs are 
equal, the product of the two type C symmetries commutes with the hamiltonians and this 
system is not minimal. More generally, considering more combinations of type P and type 
C symmetries does not lead to new minimal classes. 



B. List of classes. 



In this sub-section we present the detailed structure of the resulting classes of hamil- 
tonians. The type P and C symmetries impose the following relations on the generalized 
potentials and masses: 

Psym. : 1 A z = A zl , 1 A- = A- 1 , 7 V± + \4 7 = (2.14) 
C sym. : aA T z +A z a = $ , aAl + A-a = , a \/J = ±e c V± cr (2.15) 

They possess a simple interpretation as they indicate that A z and A- belong to an orthog- 
onal or symplectic Lie algebra depending whether cr is symmetric or antisymmetric. The 
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compatibility relations (2.8), a = ±7<7 7 T , ensures that the constraints (2.14) and (2.15) 
may be imposed simultaneously. When imposing both a type P and type C symmetry, 
one generates all relations obtained by successive applications of these symmetries; i.e. all 
relations associated to elements of the group generated by the type P and C symmetries 
are imposed. As a consequence, there could be different presentations of the same class 
depending which generators of this group one selects. For example, given a type P and a 
type C symmetry with a sign e c , their product is again a type C symmetry but with an 
opposite sign — e c . (See eq. (2.1).) 

Solutions of the constraints (2.14,2.15) for the set of compatible 7 and a give the following 
minimal classes: 

class : A-, A z G gl(n); V± G gl(n). (2.16) 

class 1 : A z G gl(n); V± = 0. (2.17) 

class 2 : A z = diag(a + , a_), a± G gl(n); 

V±= {w ± "o)' U± ' w ±^9l(n). (2.18) 

class 3 ec : A z = ^ = — A T Z G so(n), a = —a T , b = —c T ,d = —d T ; 

V-e c = ~VT tc G sofa); V ec = Vj c G gl(n) \ so(n). (2.19) 

class 4 6c : A z — ^ ^) = ~ a y^z°y e sp(2n), a = —d T , b = b T ,c = c T ; 

V- €c = -o y V\a y G sp(2n); V €c = a y V € T c a y G gl(2n) \ sp(2n). (2.20) 

class 5 : A z = -A T Z G so{n); V± = 0. (2.21) 

class 6 : A z = -a y A T z a y G sp(2n); V± = 0. (2.22) 

class 7 : A z = diag(a + ,a_), a± = — a± G so(n); 

class 8 : A z = diag(a + ,a_), a± = —a y a±a y G sp(2n); 

V± = ( w ° ± W Q ) ' ^ = ±fT ^±^" (2-24) 
class 9 ec : A z = diag(a, — a T ), a G gl(n); 

V±e c = ( Y ) ' ^ = T ^ ^ = (2 ' 25) 

(The labels 1-9 refer to the cases 1-9 listed in the previous subsection.) The hermiticity con- 
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straints A? — A\, V± — V±, v± — w±, are implicit in this list. (The hermiticity constraint is 
not made explicit for the purpose of describing the non-hermitian classes in the next section.) 
The index e c refers to one of the two possible values e c = ±; and the absence of such index 
means that this value is irrelevant. The first class corresponds to generic Dirac hamiltonians 
with no constraints imposed. There are a few degeneracies when solving eqs.(2.14,2.15). As 
expected, cases 9 and 9', eq.(2.13), yield the same solutions; also realizations of the cases 5 
and 6, eqs. (2.9,2. 10) are independent of the choice of the sign e c . Realisations of the cases 
7± are also equivalent because they correspond to different presentations of the same class. 
Indeed, consider eq.(2.11) in case 7. The product of its type P and type C defining symme- 
tries produces a new type C symmetry with opposite sign and with a' — a z <g> 1. It is gauge 
equivalent to the original symmetries as a' ~ u-j a' uj = 1 2 <8) 1 with 117 = diag(l, i)<g>l. Hence 
case 7+ and 7_ are gauge equivalent, and we give the two presentations in the above list. 
The corresponding realizations are related by similarity transformations H — > H' = U HU^ 
with U = diag(-u 7 , M7), so that a'± = ±a± and v'± = —iv±, w'± = iw±. A similar argument 
applies to the cases 8± showing again that they are equivalent presentations of the same 
class. 

One of the origins of the distinction between the above classification and the classification 
of random matrices arises from a difference in the notion of equivalent classes of compatible 
symmetries. For random matrices, some of the cases with 7 = lor7 = <r z ®l are considered 
as equivalent because they correspond to the same operator P up to re-shuffling of the lines 
and columns, while in the present classification they yield different classes because we impose 
the 2x2 block structure (1.1) to the hamiltonians 2 . Thus the classification of random Dirac 
operators is a bit finer than the one of random matrices, as summarized in Table 1. 

C. Symmetry groups, disorder measures, and super-symmetric effective actions. 

Each class is stable under the action of a symmetry group, whose elements act on the 
hamiltonians by conjugation such that their form imposed by eqs. (2. 16-2.25) is preserved. 
Elements G of the symmetry groups satisfy 3 : 

G 1 G- 1 = 1 , GaG T = a (2.26) 

The list of these groups is given in Table 1. For classes 2, 7, 8, in which the symmetry group 
is a product of two subgroups, the embedding is diagonal with G = diag((? + , g_) where g± 



2 The classification of random matrices may be reread from the previous classification by consid- 
ering that 7 and a implement directly the discrete type P and C symmetries. Only the cases with 
7 = <? z ® 1 is then relevant since 7 = 1 is trivial. As it should be, we are left with ten classes 
0, 2, 3±, 4±, 7, 8, 9-|-. However, in Table 1, the random matrix classes refer to those defined by 
C and P and not by 7 and a. See the appendix. 

3 One may extend slightly the symmetry group by discrete groups, made of Z 2 factors, by allowing 
signs in eq.(2.26), G7G- 1 = ±7, G a G T = ±a. 
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belong either to U(n), 0(n) or Sp(2n). In class 9 the embedding is G = diag(g,g T ) with 
9 e U(n). 

The symmetry group may be used to specify the disorder measures in each class, which we 
assume to be Gaussian, with zero mean, and local. The measures are then fixed by requiring 
them to be invariant under the symmetry group. The list of all quadratic invariants for each 
class is the following: 



class 
class 1 
class 2 
class 3 e 
class 4 e 
class 5 
class 6 
class 7 
class 8 
class 9 f 



tr(A z A z ), tr(A z )tr(A z ), tr(\/ ± 2 ), tr(V±) 2 , tr(V±V T ), tr(V±)tr(V T )- 

tr(A z A z ), tr(A z )tr(A z ); 

tr(a±a±), tr(a±)tr(a±), tr(v±v±), tr{v±v^)\ 

tr(A z A-), tr(Vl c ), tr(Vj 2 ; 

tr(A z A-), tr(Vl c ), tr(Vj 2 ; 

tr(A z A-); 

tr{A z A-)- 

tr(a±a±), tr(?4u±); 
tr(a±a±), tr(v ] ±v±)\ 
tr(aa), tr(a)tr(a), tr(v±v±), tr(v±v T ). 



To preserve rotation invariance we only list the invariants which couple A z to A- and V± 
to itself or to V T . Couplings between V± and V T break the symmetry under reflection 
x -> x, y -> -y 

These ensembles may be analyzed using the supersymmetric method [28] . For each class, 
this leads to an effective field theory description with the number of coupling constants equal 
to the number of invariants. These coupling constants, which measure the strength of the 
disorder, parameterize perturbations of the free field theory valid in the absence of disor- 
der. In two dimensions, all of these effective field theories can be formulated as left-right 
current-current perturbations, where the couplings are marginal [20]. This means that in the 
effective field theory the coupling constants are dual to operators of scaling dimension two. 
The discrete symmetry defining the classes plus the global invariance under the symmetry 
group should ensure that for each class the effective field theory is perturbatively renormal- 
izable, i.e. no additional marginal operators beyond those dual to the coupling constants are 
generated by the renormalization procedure. This aspect can be studied using the all-orders 
P function proposed in [29], as was done for class and for class 4_ at n = 1 [30]. 

We can easily describe the global Lie superalgebra symmetry of the effective field theories. 
The unperturbed conformal field theory has an osp(2N\2N) current algebra symmetry at 
level 1 where iV is the number of fermions, i.e. iV = n for classes 0, 1, 3 and 5, N = 2n for 
classes 2, 4, 6, 7 and 9, and N = 4n for class 8. In the supersymmetric effective theory the 
global supersymmetry is smaller, and corresponds to the Lie superalgebraic extension Q of 
the symmetry groups listed in Table 1. The bosonic group U (n) is extended to gl(n\n), 0(n) 
to osp(n\n) and sp{2n) to osp(2n\2n), so that Q = gl{n\n) for classes 0,1,9, Q = osp{n\n) 
for classes 3,5, Q = osp(2n\2n) for classes 4,6 and a tensor product of these supergroups 
for classes 2,7,8. 
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III. CLASSIFICATION OF NON-HERMITIAN DIRAC HAMILTON! ANS. 



The classification of non-hermitian Dirac hamiltonians we present is based on imple- 
menting discrete symmetries of type P or C, eqs.(1.2,1.3), and of type Q and K defined in 
eqs.(1.5,1.4). Since this parallels closely what we have done for hermitian hamiltonians we 
only sketch the main features of the classification. The order two constraints on type Q and 
K symmetry are: 

& = 1, e = 1 ; KK t = 1, k t = ±k 

We consider hamiltonians up to unitary changes of basis, H — > UHW, with U = diag(w, u), 
which act on Q and K as: 

^^>u^v) ; k-^uku t (3.1) 

As before we define minimal classes as those whose hamiltonians do not commute with a 
fixed matrix preserving their Dirac structure. 

Imposing the type Q or K symmetries amounts to imposing some reality conditions on 
the hamiltonians, i.e. some reality properties of A z , A- and V±: 

Qsym.: £Al = A z £ , ZV± = e q V±Z (3.2) 
Ksym. : kA^+A z k = , «VJ = ±e fc V r ±« (3.3) 

Redefining H — > iH modifies the signs e q and e& in eqs.(1.5,1.4), however this redefinition 
ruins the Dirac structure (1.1) and we shall thus not allow it. 

The classification of non-hermitian hamiltonians may be translated into detailed proper- 
ties of the hermitian hamiltonians Ti. obtained by doubling the Hilbert spaces on which the 
Dirac hamiltonians H are acting: 

«=(£ o) < 3 ' 4 » 

These doubled hamiltonians are always chiral as they anticommute with T 5 = diag(l, — 1). 
Any similarity transformation H — > UHU' 1 is mapped into H — > UHU^ with U = 
diag(C/, U^ 1 ). Demanding that these transformations also act by similarity on Ti, imposes 
U to be unitary. When no discrete symmetries are imposed, the doubled hamiltonians Ti 
are always elements of class 2, which is embedded in the chiral GUE class. Indeed, up to 
re-shuffling of lines and columns, they may be presented as: 



H~H d 



( Vl + Vl -idz + Al \ 

V + + V_ -idz + A- 

-id z + A z V+-V- 

V -id g + Al Vl-Vl ) 



(3.5) 



The dictionary is thus a + = A z , a_ = A% and 2v± = (Vj ± V+) + (Vl =F VL). 

On H : both type P and Q symmetries act as chiral transformations, H — > —VHP -1 

with P = diag(P,P) and H -> e 9 Q7iQ _1 with Q = (g Thus > W may be block 
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diagonalized if H is Q or P symmetric. Indeed, if H is Q symmetric with e q = +1 then Q 
and TL may be simultaneously diagonalized since they commute. Similarly, if H is P or Q 
symmetric with e q = — 1, then TL commutes with the product T 5 V or with T 5 Q. 

Type C and X symmetries both act as particle-hole symmetries relating TL to its trans- 
posed TL T . The classification of the hamiltonians TL is then very simple as it follows from 
that of hermitian Dirac operators. For the operators TL to be minimal, only a type C or 
a type K symmetry can be imposed. Gauge equivalences (2.2,3.1) leave only a = 1 or 
a = i<j y <S> 1 and k — 1 or k — ia y <g> 1 as possible choices. We then have the correspondence: 

a = 1 => TL E class 7 (3.6) 
a = io y ® 1 =>- TL G class 8 
k = 1 or k = io y (g> 1 =^ 7i G class 9± 

Though we have translated the classification of non-hermitian H into the doubled her- 
mitian TL, the spectrum of TL and H may differ^ significantly. To illustrate these potential 
differences, consider the transformation H — > H = —iujHuj where 117 = diag(l,i) ® 1 is 
the matrix we introduced in Section II. This transformation preserved the Dirac structure 
of the hamiltonians but not their reality conditions. It leaves invariant A z and A- but not 
the potentials since V± = —iV T . Hence H and H should not belong to the same non her- 
mitian class — and they do not have the same spectra. On the contrary, for the doubled 
hamiltonians this transformation is lifted to TL — > H = UTtlfi with U = diag(— ill?, u\), so 
that TL and TL are unitarily equivalent. They have identical spectra and belong to the same 
class. 

We thus present a classification of the Dirac operators H and not simply of the doubled 
one. As for hermitian Dirac operators, it is the group generated by compatible discrete 
symmetries which is meaningful. There could be different but equivalent presentations of 
the same group as not all of these symmetries are independent. Indeed, the product of a 
type P symmetry with a C, K or Q symmetry is again a C, K or Q symmetry. (See for 
example eq. (2.1).) Also, the symmetries of type Q, C or K are linked as the product of 
any of two of them gives a symmetry of the third type. 

Let us first impose only one type of symmetry. As in Section II, up to gauge equivalence 
(2.2,3.1), the solutions are: 

(7 = 1) ; (<7 = ik ; (£ = ik ; = ; 

( 7 = a z ® 1) ; (a = ia y <g> l) £c ; (f = a z <g> l) £q ; (k = %a v <g> \) tk ; 

Here, each column refers to one of the possible types of symmetry. Here and below, we 
indicate as indices the values of e c , e q or which matter. Thus the above list corresponds 
to 14 distinct classes. 

Let us now impose two kinds of symmetry. First, we may require simultaneously type 
P and C symmetries. This leads to the list of 6 classes, from class 5 to class 9 of Section II 
without any reality conditions. 

Next we consider a P and a K symmetry. The commutativity condition for type P 
and K symmetries reads k = ±7~ 1 /t7*. This is solved the same way as o = ±7 _1 o"7 T 
in previous Section. Thus the list of compatible type P and K symmetries is parallel to 
the list of compatible type P and C symmetries, only a is replaced by k. Their explicit 
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realizations are given in eqs. (2.21-2.25) but with v± replaced by w± and w± replaced by v±. 
This corresponds to six classes. 

Similarly compatibility between type P and Q symmetries requires 7^ = ±£ _1 7 £. Solving 
this constraint leads to the following compatible type P and Q symmetries: 

(7 = 1, f = 1) ; (7 = 1 <g> 1, f = a z <g> 1) ; (7 = a z <g> 1, f = <r x <g> 1)^; 

(7 = <7* <g> 1, £ = 1 ® l)e q = (7 = ^ ® 1, £ = ^ ® l)- e , S 

In the second line, we have mentioned an equivalence between two solutions of the commu- 
tativity constraint. Indeed, £ of the second solution in this line is the product of 7 and £ of 
the first solution, so the groups generated by these solutions are identical. 

We may also impose together a type Q with a type C symmetry. Since their product 
is a symmetry of type K with = e q e c , we are actually imposing simultaneously three 
compatible symmetries of different types. Any two of them generate the third. The condition 
for the type Q and C symmetries to commute is £ T = ±er £ _1 a, for type C and K this 
condition reads k t <j~ x ku* = ±1. Up to gauge equivalences, the set of compatible type Q 
and C symmetries is then: 

(f = 1, a = l) eq , ec ] (f = a z <g> 1, a = 1 <g> 1)^; (f = <r z <g) 1, <7 = ia y <g> 1)^; 

(f = 1 <g> 1, (J = Wy <g> l) eq , ec ; (f = <7 Z <8) 1, (J = 1 <8 iVy)e q ,e c ; (£ = 0" 2 (g) 1, (7 = O x <8 l) eq , £c 

Finally, we may impose simultaneously a type P symmetry together with two among the 
three types Q, C and K of symmetries. As before it is sufficient to consider only a P,C and 
Q symmetry. The solutions of the commutativity requirements are then: 



(7 = 


18 


'1, t = 


1 8 


1, <r = 


= 1 <8 


1 or itjy ® 1); 




(7 = 


18 


»i, i = 


(T 2 C 


3 1, <7 


= U 


5 1 , cr = 1 <8> icry , ^0"^ ®\ oy a x (i 


5 1); 


(7 = 


a z ( 


8>i, i-- 


= lc 


g) 1, (7 


= U 


5 1 or 1 ® i(Jy)e q ; 




(7 = 


a z ( 


8)1, £ = 


= lc 


g) 1, a 


= o"x 


® i) e , )ec ; 




(7 = 


a z ( 


8U = 




8> 1,(7 


= l€ 


) 1, 1 ® icT^, cr s ® 1, or (j x ® io y £ 


3l) e , 



Here, each choice of a corresponds to a different class. 

Considering more combinations of the four different kinds of symmetries would not lead 
to new minimal classes, because in such case we would always be able to construct some 
matrix commuting with the hamiltonians and preserving their Dirac structure. 

For each set of compatible symmetries, one has to choose the signs e q , e c and e& to specify 
the classes. These signs are used to index the solutions in the above lists. The absence of 
one of this index means that the corresponding solution is independent of that index. The 
grand total is 87 classes. 

It is straightforward to determine the form of H for each of the above symmetry classes, 
however there is little motivation to list the details here. Let us just describe a simple 
example, corresponding to imposing a symmetry of type K, relating H to its complex 
conjugate, with k — 1 and = +1. Then, relations (3.3) yield A- = —A* z and V± = ±V^, 
such that V + is real and V_ imaginary. As a consequence, the Dirac hamiltonian may be 
written as: 

( M -ids-Al\ 
\-id z + A z M* J 
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with M — V + + V-. This class, studied e.g. in ref. [31], is closely related to the random XY 
model. The doubled hamiltonian belongs to class 9 . 

Having performed the above classification we can easily specialize it to random non- 
hermitian matrices with no Dirac structure. As for random hermitian matrices (see the 
appendix for a summary), the above classes with 7 = 1 are trivial and should be thrown 
away. The choice of the sign e q and e& is also irrelevant since it can be absorbed into H — > iH 
which is now allowed since no Dirac structure is imposed. Altogether this gives 43 classes 
which will be described in greater detail in [32]. 

IV. APPENDIX. 

For completeness — and for explaining Table 1 — we recall the definition of random 
hermitian matrix ensembles [15]. We denote by small letters quantities referring to random 
matrices. Let h = h) be a hermitian matrices and p and c be the operators implementing 
the discrete symmetries as in eqs.(2.14,2.15): h — > —php -1 , h — ■» e c c/i T c _1 . For each random 
matrix class, the defining relations for p and c are summarized in Table 2. 



Random matrix 
classes 


Discrete sym. 
Relations 


A 


h = W 


A, 


c T = c, e c = + 


A n 


c T = -c, e c = + 


A/// 


p 2 = l 


C 


c 1 = -c, e c = - 


D 


c 1 = c, e c = - 


D 7 


p 2 = 1, c 1 = c, e c = ±, pep 1 = c 


C77 


p 2 = 1, c 1 = — c, e c = ±, pep 1 = c 


Cj 


p 2 = 1, c 1 = ±c, e c = ±, pep 1 = —c 


D77/ 


p 2 = 1, c T = ±c, e c = =F, pcp T = -c 



In each of the last four lines of Table 2, one may equivalently choose either the upper 
or the lower signs, since this simply corresponds to choosing two equivalent presentations of 
the same class. 

To compare with the classification of Dirac fermions, it is useful to notice that in the 
latter case the operators P and C may be written as: 

P = r z <8> 7; C = ir y <S> o for e c = +; C = r x <8> a for e c = — . 
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